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We study to what extent the detrimental impact of dissipation on quantum properties can be 
compensated by suitable coherent dynamics. To this end, we develop a general method to determine 
the control Hamiltonian that optimally counteracts a given dissipation mechanism, in order to 
sustain the desired property, and apply it to two exemplary target properties: the coherence of a 
decaying two-level system, and the entanglement of two qubits in the presence of local dissipation. 
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Genuine quantum features such as entanglement or co- 
herence are resources as precious as fragile, and their 
uncovering usually requires strong efforts in isolating 
and controlling quantum systems. Without thorough 
measures, decoherence efficiently shields the quantum 
world from our access and hides it behind its classical 
guise. While there has been unprecedented progress in 
the quantum control of various model systems, e.g. ions 
PQ , quantum dots [5] , or cold atoms [3] , it is not possible 
to completely decouple these systems from their environ- 
ment and thus to fully suppress the detrimental effect of 
decoherence. Standard optimal control techniques there- 
fore focus on accessing quantum features in the transient 
regime, and the exploration and exploitation of quantum 
properties is consequently confined to a finite, generically 
short time window. 

There are however ways to keep the window to the 
quantum world enduringly open, e.g. by encoding quan- 
tum features in topological properties of a system [3] , or 
by engineering a dominant environment that drives the 
system into a highly nonclassical stationary state [6] . 
While these approaches in principle permit to prepare ar- 
bitrary nonclassical quantum states, they require in gen- 
eral an exceedingly large overhead of resources. 

In this paper we therefore ask to what extent already 
standard Hamiltonian control can enduringly counteract 
the detrimental effect of decoherence. Explicitly, we seek 
Hamiltonians that optimally uphold, on asymptotic times 
scales, a given control objective (e.g., coherence, entan- 
glement, or fidelity w.r.t. a target state) in the presence 
of dissipation. Such an asymptotic time behavior can 
be meaningfully formulated for static and, more gener- 
ally, periodically time-dependent Hamiltonians. In the 
latter case, the asymptotic dynamics are periodic cycles 
in state space, which reduce to stationary states in the 
static case. 

With respect to our goal to single out the optimal 
among all conceivable control Hamiltonians, it is not ad- 
visable to directly scan the space of Hamiltonians, as the 
latter cannot efficiently be parametrized. We therefore 
approach the problem from a different perspective and 
determine the optimal stationary state or asymptotic cy- 



cle directly, i.e. independently of the Hamiltonian. The 
crucial insight behind is that physically admissible trajec- 
tories in state space are strongly constrained by the dis- 
sipative part of the dynamics. It thus turns out that one 
can characterize all possible stationary states or asymp- 
totic cycles from the dissipative dynamics alone. 

While our approach can be applied to the optimization 
of arbitrary control objectives, we demonstrate its viabil- 
ity with two physically relevant examples: the coherence 
between ground and excited state of a decaying two-level 
system, and entanglement of two qubits in the presence 
of local dissipation. In both cases, we find that the opti- 
mal Hamiltonians uphold substantial, but limited values 
of the respective control objective. 

Static control Hamiltonians. Throughout the paper, 
we consider an open quantum system that evolves under 
a Lindblad master equation [7], 

p(t) = i[p(t),H(t)}+V(p(t)), (1) 

with a dissipator V(p) = Y.klk[ L kpL\ - \{L,\h k , p} + ] 
composed of Lindblad operators L k and rates jk- For an 
arbitrary but fixed dissipator T)(p), our goal is to opti- 
mize the stationary state p ss (for static H) or asymptotic 
cycle p ac (t) (for periodic H{t)) of ([I]) with respect to an 
arbitrary objective function O(p). 

It is instructive to investigate static control Hamiltoni- 
ans first. A direct optimization in the set of all conceiv- 
able Hamiltonians H involves the stationarity condition, 

= i[ Pss ,H] +V{ Pss ), (2) 

in order to infer the stationary state p ss associated with 
a particular choice of H. Its inversion typically requires 
numerical means and must be repeated for each sam- 
ple Hamiltonian, rendering this approach impractical al- 
ready in low-dimensional quantum systems. Therefore, 
we develop a different strategy here. Instead of taking 
Hamiltonians as our starting point, we base the opti- 
mization on the set of stabilizable states S [5]: 

S = {p: 3H s.t. = i[p, H] + V(p)}. (3) 

By definition, it comprises all quantum states that be- 
come stationary under a suitable Hamiltonian. As shown 



below, this set can be characterized independently of the 
specific Hamiltonian. Optimization of an objective func- 
tion O(p) can then be done in S directly. 

To derive this Hamiltonian-independent characteriza- 
tion of S, we exploit that the coherent dynamics induced 
by H and the dissipative dynamics induced by V(p) must 
compensate each other for a stationary state. Since the 
coherent part of the master equation ([!]), i[p,H], neces- 
sarily leaves the spectrum of p invariant, this must also 
hold for the dissipator T>(p) at a stationary state. In 
other words, T>(p) does not modify the purity p — Tr[p 2 ] 
or any higher moment Tr[// 1 ] with n > 2. This implies 

9tTr[p n ]\ H=0 S nTi[p n - 1 V(p)] = 0, if p is in S. Since 
the moments of a d-dimcnsional quantum state are inde- 
pendent only up to n = d, this leads to d — 1 necessary 
conditions for p to be stabilizable: 



(a) 7+/7_ = 7d/7- = (b) 7+/7- = 7d/7- = V 4 



pes 



Vne{2,...,d}: Ti[p n - 1 V(p)}=0. (4) 



Denoting by S n the set of states that fulfill criterion Q 
for fixed n, we have S C f~) 5 n . For states with non- 
degenerate eigenvalues, criterion Q is also sufficient (see 
supplementary material) , so that we have S — |~) n S n in 
this case. Note that this hierarchical characterization of 
S does not require reference to the stabilizing Hamilto- 
nian H - in contrast to the original definition ([3]) . Given 
a stabilizable state p £ S, however, it is straightfor- 
ward to derive the corresponding H from the station- 
ary condition pi), based on the spectral decomposition 
P = T ia ^a\a) (a\: 



H = 



E 

a,/3:A Q #A,3 
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To demonstrate the viability of our method, we first 
discuss the case of a single qubit. There, one can find 
an intuitive geometric representation of the set S. Since 
d = 2 in this case, criterion Q imposes merely a single 
constraint (n = 2). In terms of the Bloch vector, f — 
Tr[pa], this constraint defines a quadric hypersurface in 
the Bloch ball: 



f G S 



(Dr + d\ 



d) = 0. 



(6) 



Here, the 3x3 matrix (-D)jj = Tr[<jjD(cry)] and the con- 
stant vector (d)i = Tr[<7jX>(l)] characterize the dissipator 
in Bloch notation. According to ([6]), a state r is stabi- 
lizable, if (and only if [20]) the dissipative flux Df + d is 
orthogonal to r, i.e., if it has no radial component. 

To be specific, we consider a qubit that is exposed to 
the three most common incoherent processes: decay of 
the excited state at rate 7_ , absorption from the ground 
state at rate 7+, and dephasing between ground and ex- 
cited state at rate 7^. In this case, S is the surface of a 
spheroid [8], with the polar axis pointing in z-direction, 
cf. Fig. [l] The polar and equatorial diameter depend on 




FIG. 1: (color online). Bloch ball representation of the set 
S of stabilizable states (surface of the orange (dark shaded) 
spheroid) for a single qubit that is subject to spontaneous 
decay, excitation, and dephasing. The respective rates 7_, 
7+, and 7^ determine the shape of S. If only spontaneous 
decay is present (left panel), the stationary state r» that is 
optimal with respect to the coherence C (i.e., the distance 
to z-axis) reaches C = l/v2. With finite 7+ and 7^ (right 
panel), the ellipsoid is contracted and the optimal value of C 
is reduced. 



the incoherent rates. The optimal stationary state with 
respect to an arbitrary objective function 0(r) is now 
conveniently determined by maximization over the sur- 
face of this spheroid. E.g., one may consider the coher- 
ence C — 2| (0|p|l) I between ground and excited state of 
the qubit. In Bloch notation, this objective can geomet- 
rically be interpreted as the distance to the z-axis, see 
Fig. [T] Hence, the optimal coherence equals the equa- 
torial semi-axis of the spheroid, which is r_/fi, where 
n = v/r + (r + /2 + 7 d ) and r± = 7_ ± 7+ According 
to Eq. ([5]), the Hamiltonian that stabilizes the optimal 
stationary state r* is then H* = —^cx y . It can be re- 
alized, e.g., by resonantly driving the qubit with a Rabi 
frequency f2 [S]- 

As an example for a higher-dimensional quantum sys- 
tem, we consider two qubits (d = 4). The set of stabi- 
lizable states S lies in the intersection of three hypersur- 
faces 1S2, S3, and £4 in this case. Similarly to a single 
qubit, one can represent the two-qubit state p by a gen- 
eralized, 15-dimensional Bloch vector r [TU]. The lowest 
order constraint (n — 2) can then be cast into the same 
form as ([6]), defining again a quadric surface for r. The 
higher order constraints for n = 3,4, however, lead to 
polynomial expressions of third and fourth degree in r. 
Therefore, instead of determining the optimal state in S 
directly, it is favorable to determine the optimal state F* 
in S2 first, and then to verify that r* lies also in S. If 
so, it must be the optimum in S, since 5 C If not, 
the procedure provides at least an upper bound for the 
optimal value in S. 

Relevant target properties for two qubits address, e.g., 
their entanglement. To this end, it suffices to solely con- 
sider spontaneous decay at rate 7_, because additional 
processes like absorption or dephasing merely contract 
the quadric surface S2 towards the totally mixed state 
r = 0, cf. Fig. [lj This corresponds to depolarizing every 
state in S by a factor x (i.e., r — ¥ xr with < x < 1), 
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implying that the optimal entanglement in S is always 
reduced in the presence of these additional processes. 

As a specific entanglement objective, we study the fi- 
delity F{p) = ($+|p|W+) with respect to the maximally 
entangled Bell state |*+) = (|01) + \\Q))/y/2. It is re- 
lated, e.g., to the fidelity of a teleportation protocol via 
the state p [11 . In this case, the optimization in 52 can 
be carried out analytically, leading to the optimal sta- 
tionary state 

p* = l/2|Q0)<00| + l/2|¥ + ><¥ + |, (7) 

with fidelity F{p*) = 1/2. The Hamiltonian H(a,j3) = 
l®(aa z + /3a x ) + (aa z + /3a x )(g>l-2a(a + (g>a_+a^<E)a + ) 
stabilizes p* in the limit of f3/j_ — > oo and a/ ft — > oo 
[21j . This optimal Hamiltonian is readily realized in var- 
ious experimental setups; in particular, the interaction 
term corresponds to an excitation hopping mechanism 
that describes, e.g., dipole-dipole interactions between 
excitons [TH [T3] or Rydberg atoms [21 [TS] . 

Another relevant two-qubit objective is the entangle- 
ment measure concurrence £{p) |16j . In contrast to the 
fidelity, it does not favor a specific state, but it assigns the 
full value £ = 1 to all maximally entangled states. Since 
£(p) is not a linear function of p, however, its optimiza- 
tion can no longer be treated analytically. Numerical op- 
timization in £2 reveals that the optimal state coincides 
with the fidelity- optimal state p* of Eq. ([7]), which car- 
ries a concurrence of £ = 1/2. This finding is plausible, 
since a high Bell state fidelity typically goes along with 
strong entanglement. With respect to the concurrence, 
however, p* is not the unique optimum; e.g., ^ with ty + 
replaced by ^_ is a stabilizable state with concurrence 
1/2, as well. 

These results clearly indicate the possibilities and lim- 
itations of coherent control of open systems: on the one 
hand we have proven that it is not possible to get be- 
yond the fifty-fifty mixture Q of the deexcited state and 
a maximally entangled Bell state; on the other hand, 
£ = 1/2 is still significantly higher than the average con- 
currence £ = 0.18 of typical stationary states, as we have 
verified by a statistical sampling of Hamiltonians. 

Periodic control Hamiltonians. So far, we developed 
a method to determine the static Hamiltonian H that 
upholds the optimal amount of a given objective O in 
the stationary state. The derived concepts can however 
be extended to the envisaged, substantially more general 
case of periodic control Hamiltonians H(t) = H(t + T) 
with period T, Since this comprises static Hamiltonians 
as a special case, one generally expects that this addi- 
tional freedom in control should improve the optimiza- 
tion results. 

Specifically, our aim is to determine periodic control 
Hamiltonians H(t) which optimize the objective function 
of the resulting asymptotic cycle p ac (t) on time-average, 
i.e., which maximize O = |/ 0(p ac (t))dt. To this end, 



we optimize O in the set A of stabilizable cycles, com- 
prising all periodic trajectories p(t) = p(t + T) (with 
arbitrary period T) for which a periodic Hamiltonian 
H(t) — H(t + T) exists such that the trajectory becomes 
the solution of the master equation ([I]). Criterion Q is 
then generalized to 

Vt Vn : TT[p n - l {t)V{p{t))] = -d t Ti[p(t) n ], (8) 

n 

which must hold for any {p(t)} € A. Eq. {[sj) reflects the 
fact that only the dissipative term T>(p) can change the 
spectrum of p and thus its spectral moments Tr[p n ], but 
not the Hamiltonian part of the dynamics. As before, 
criterion Q is also sufficient for {p{t)} G A, if p(t) has 
non-degenerate eigenvalues for all t. The Hamiltonian 
H(t) that stabilizes a given cycle in A is found analo- 
gously to ([5]). 

E.g., for the purity p (i.e., n = 2), criterion ([8| im- 
plies that a consistent cycle must equally probe regions 
of the state space where the purity flux f(p) = Tr[pT>(p)] 
of the dissipator is positive and regions where it is neg- 
ative. The regions of positive and negative purity flux 
are separated by the hyperplane of vanishing flux, which 
coincides with the set 52 that we introduced above when 
studying static Hamiltonians. Hence, any cycle must in- 
tersect 1S2 an even number of times, but at least twice. 
Moreover, no cycle can explore regions where the purity 
is larger than the maximal purity in £2 (see supplemen- 
tary material). These statements hold likewise for the 
higher moments Tr[p™]. 

It is again instructive to consider the case of a single 
qubit first. A cycle {p(t) : t € [0,T)} is then represented 
by a closed trajectory {f t } in the Bloch ball, and stabiliz- 
able cycles are characterized by criterion (|Sj) with n = 2. 
In Bloch notation, this criterion reads 

rv (Dr t + d) = \d t \r t \ 2 (9) 

= /( ? t) =P{r t ) 

and reflects the fact that the time evolution of the purity 
p = Tr[p 2 ], which corresponds to (|r| 2 + l)/2 in Bloch 
notation, is exclusively governed by the radial part /(r) 
of the dissipator. 

While an optimization over all stabilizable cycles is 
certainly unfeasible, the problem can be reduced to the 
tractable class of two-point cycles (TPC), based on the 
following argument: Any stabilizable cycle must undergo 
subsequent stages of strictly monotonic purity gain and 
loss, in order to be closed. Without loss of generality, we 
consider a cycle that consists of two stages, intersecting 
S2 twice; the analysis of general cycles reduces to this 
case. To each point r^ in the purity-increasing stage 
(/(^p ) > 0)' one can assign a point of equal purity p 
in the purity-decreasing stage (f(f~) < 0). Hence, the 
cycle can by parametrized by the purity p, as sketched 
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FIG. 2: (color online). Illustration of a general asymptotic 
cycle {f t } (dotted loop), and a two point cycle {(r£ , f^)} (red 
boxes). White rings connect points of equal purity p. Orange 
(dark gray) color indicates the region of positive purity flux. 
It is separated from the region of negative flux (light gray) 
by the hypersurface 1S2 (solid black line) which comprises all 
states that can be stabilized by a static Hamiltonian. 



in Fig. [2] Denoting by po (pi) the minimal (maximal) 
purity of the cycle and using ^ , we find that the average 
objective value O of an arbitrary cycle can be bound from 
above by a TPC, 



O 



fPi ( Q(f* 

ho \\f(r + 



) , O(r-) 



)l 1/(^)1 



rPi ( l_ 

ho \\f(4 



)l ' l/(r>)l 



dp 

— < max0 TPC (rj",rt"). 
dp p 

(10) 

Here, the average objective value Otpc of a TPC is given 
by 

TPC(p,pj 1/(^)1 + 1/(^)1 ■ 1 } 

Inequality ( p~0| ) is based on the estimate 
J g(x)dx/ w(x)dx < max x g(x)/w{x), which holds 
for any positive function w{x) and any g(x). The cycle 
that achieves 0tpc(^ consists of only two points 
r+ and fp of equal purity p, lying on different sides of 
the hyperplane iS>2- To see this, one realizes that the 
(infinitesimal) purity 8p is lost while residing for an 
(infinitesimal) dwell time St" at fp~ . To close the cycle, 
this purity 5p must be regained during the dwell time 
5t + at r£ . Since the ratio of these dwell times is inverse 
to the ratio of the respective purity fluxes, 



St+ j|} 8p/\f(r+)\ = |/(r-)| 
St- 6p/\f(ft)\ \f(4)l 



(12) 



the time-average O of the TPC is given by Eq. ( 11 ). The 



TPC rapidly jumps back and forth between fj and fp 
via purity-preserving, unitary "kicks" , which can be gen- 
erated by suitable, ideally (5-shaped pulses in the Hamil- 



tonian H(t). In conclusion, inequality (10 1 reflects the 



important result that the search for the optimal asymp- 
totic cycle can be restricted to simple TPCs, simplifying 



tremendously the original optimization over all closed 
trajectories that are consistent with criterion 

To exemplify this reduction, we consider again the co- 
herence C of a single qubit undergoing the same incoher- 
ent processes as before (characterized by the rates 7_, 
7 + , and j z ). Due to symmetry reasons, the most gen- 
eral TPC is then parametrized by two azimuthal angles 
and the common purity p, and it is constrained by the 
condition /(fp") > 0, f(f~) < 0. The optimization of 
^TPcirp ,f~) over all TPCs is then readily performed 
numerically. One finds that for any combination of the in- 
coherent rates, the optimal TPC {f^,f^} degenerates to 
a single point, namely, to the optimal stationary state r* 
for static Hamiltonians that is marked in Fig. [TJ Remark- 
ably, this implies that no periodic Hamiltonian H(t), no 
matter how sophisticated, can beat the optimal static 
Hamiltonian H* . This is, however, a peculiarity of our 
specific choice of objective and dissipator and does not 
hold in general. 

The same strategy can equally be applied beyond a 
single qubit. Similar to the static case, one obtains an 
upper bound for the optimum in A by focusing on the 
lowest-order set Ai alone, since A C Ai- This way, 
one can again restrict the investigation to TPCs, render- 
ing a numerical optimization feasible. In our two-qubit 
example, we find that the time-averaged concurrence 
£ = ij; Jq £{p(t))dt never exceeds the optimal static re- 
sult £{p*) — 1/2, if only spontaneous decay is present. 
This value decreases in the presence of a finite absorption 
rate 7+ and/or dephasing rate 7^, and it never outper- 
forms the static optimum. 

Conclusion. We developed a method to characterize 
the asymptotic states of open quantum systems in terms 
of the dynamical constraints imposed by the dissipator. 
It allows us to access the optimization of arbitrary peri- 
odically time-dependent coherent control without resort- 
ing to the system Hamiltonian. In the static case, the 
method leads to the characterization (4) of stabilizable 
states, which reflects the unitary compensability of the 
dissipator. Based on that, we showed that in the gen- 
eral periodic case, optimizations can be restricted to the 
significantly simplified class of two-point cycles. This 
way, optimization problems w.r.t. arbitrary objectives 
can be addressed that were previously prohibited by the 
vast range of conceivable asymptotic cycles. To prove 
the viability of our approach, we determined the max- 
imum asymptotic two-qubit entanglement that can be 
preserved by periodic coherent control in the presence 
of a dissipation- inducing environment. Other relevant 
applications include, e.g., optimal energy transport in 
quantum networks [IT] or the minimization of particle 
loss in Bose-Einstein condensates [18]. Altogether, our 
method not only deepens our conceptual understanding 
of the working principles in open quantum systems, but 
also opens the prospect to treat hitherto intractable op- 
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timization problems. 
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The "only-if" holds for any state apart from the totally 
mixed state r = 0, which has degenerate eigenvalues, so 
that criterion |6| is not sufficient for stabilizability. How- 
ever, lies in the closure of <S, as apparent from Fig. [l] 
We strongly conjecture that f") S n is always identical to 
the closure of S, since states with non-degenerate spec- 
trum are dense in the state space, making it plausible 
that they are dense in S as well. 
[21] Thus, strictly speaking, p* is not in S, but in its closure, 
consistent with the fact that is has degenerate eigenval- 



In the following, we prove that criterion Q is sufficient 
for pGS. Let p be a quantum state on a d-dimensional 
Hilbert space. If p has non-degenerate eigenvalues, and 
if 



Tt[p n -'V(p)] = 



(13) 



holds for n — 2, . . . ,d, then there is a Hamiltonian H 
that renders p the stationary state of 0, i.e., p € S. 

Proof: By assumption, p — J2 a ^« \oe) (a\ has non- 
degenerate eigenvalues A Q . In order to prove that p is 
stationary under 0, we show that 

(a\D(p)\P)=i(a\[H,p]\0) Va, (3 = 1, . . . , d, (14) 

with the Hamiltonian H of Eq. Since {\a}} is a com- 
plete basis, expression ( 14 ) is equivalent to the operator 
identity = i\p, H\-\-T>(p), and hence implies stationarity 
of p. 

Inserting ^ into (14), and using the notation D a p = 
(a\T>(p)\/3) , we have 



= i 



iDg'f}' 

, A Q ' — A 



{a\(\a>) (f3'\p-p\a') tf\)\P) 



An 



Xr- 



(Note that we could replace X a i ^ \p> by a' ^ f3' in the 
summation, due to the non-degeneracy of the eigenvalues 
{A a }.) The above expression is obviously true for a ^ 
(3. It remains to show that D aa = holds for all a — 
, d. To this end, we rewrite condition ( 13 1 like 



1 



= Tr[p n - 1 V( P )} = J2(K) n - 1 D o 



(15) 



By assumption, this holds for all n = 2, . . . , d. In addi- 
tion, it also holds for n — 1, since any dissipator fulfills 
Tr[2?(p)] = 0. (This can be directly verified from the 
Lindblad form of T>). Hence, one can write (151 as a 
matrix equation SDto = 0, with (DJV) an = (A a )"^, and 
(o) a = D aa . dJl is a Vandermonde matrix, for which 
det9H = ricK/sC^a — A^) is known [TS]. Since the eigen- 
values A a are non-degenerate by assumption, we have 
det9H 7^ 0, and 93? is hence invertible. The only solution 
to QJt = is therefore = 0. Thus, we have shown that 



for all a = 1, 



Supplementary Material II 

We demonstrate that no cycle can enter regions where 
the purity is larger than the maximal purity in S2- Let 
Pi be the maximal purity in the set 1S2, 



pi = max[Trp 2 ]. 



(16) 



Since 52 separates the regions of positive and negative 
purity flux f(p), f(p) has the same sign for all p with 
Tr[p 2 ] > pi- To fix the sign, we exploit the fact that 
the purity flux is non-positive for pure states, charac- 
terized by maximal purity Tr[p 2 ] = 1 - otherwise, the 
dissipative flux would increase the purity beyond 1, leav- 
ing thus the space of valid quantum states. Hence, for 
any state with purity larger than p\ , the dissipative flux 
is negative. Now, suppose that a cycle {p(t)} has purity 
Tr[p 2 (i )] > Pi at some time t . Then, at the instance 
before to , the purity must have been even larger, because 
the purity flux is negative in this region. Going back in 
time, the purity will always strictly increase, never leav- 
ing the region of negative flux, and never returning to 
its value at time to, in contradiction to the assumption 
that {p(t)} is a cycle. This proves that states with purity 
larger than p\ can never be traversed in an asymptotic 
cycle. 



